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Abstract 

We demonstrate that the integrated gravitational wave signal of Type la supernovae (SNe la) in 
the single-degenerate channel out to cosmological distances gives rise to a continuous background to 
spaceborne gravitational wave detectors, including the Big Bang Observer (BBO) and Deci-Hertz 
Interferometer Gravitational wave Observatory (DECIGO) planned missions. This gravitational 
wave background from SNe la acts as a noise background in the frequency range 0.1 - 10 Hz, which 
heretofore was thought to be relatively free from astrophysical sources apart from neutron star 
binaries, and therefore a key window in which to study primordial gravitational waves generated 
by inflation. While inflationary energy scales of > 10^^ GeV yield inflationary gravitational wave 
backgrounds in excess of our range of predicted backgrounds, for lower energy scales of ~ 10^^ 
GeV, the inflationary gravitational wave background becomes comparable to the noise background 
from SNe la. 

PACS numbers: 04.30.Db 04.40.Dg 97.60.Bw 97.60.-s 
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INTRODUCTION. 



Type la supernovae (SNe la) are powerful thermonuclear explosions of a white dwarf star 
that release a kinetic energy of ~ 10^^ ergs, and serve as standard candles for cosmology 1- 
Isj. While the precise nature of the SNe la progenitor remains a topic of active investigation, 
current theoretical models j4-6| of the detonation of a white dwarf star brought to a mass 
near the Chandrasekhar limit by a companion red giant or main sequence star are consistent 
with both the energetics and nucleosynthetic yield of the explosion, as well as its observed 
asymmetry |7[. We refer to this particular mechanism of SNe la as the single degenerate 
channel (SD)|8|, loj, as opposed to the double degenerate (DD) channel, which consists of the 
merger and subsequent detonation of two degenerate white dwarfs [lO, [ill ■ 

Two models for detonations in the SD channel have been studied extensively - the 



deflagration-to-detonation transition (DDT) 12|, ll3| , and the gravit at ionally- confined deto 



„at,o,. (GCD). Utilizing detailed three-ditnensipal (3D) hydrodynam.cal simulations of the 

GCD mechanism of the SD channel of SNe la, [14l| demonstrated that the intrinsic asymme- 
try of the explosion gives rise to a gravitational wave (GW) signature. Based on the GCD 
model prediction, a nearby individual SNe la event may possibly be detected in the Milky 
Way and nearby galaxies by third-generation GW detectors in the ~ 2 Hz frequency range. 
While the focus of the earlier paper [l^ was upon the GCD model for detonation in the SD 
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channel, both the DDT and DD models are likely to yield off-centered detonations 
Recent work has demonstrated that the explosion energies and nucleosynthetic yields are 
similar in off-centered detonations in which only a small fraction of the total energy release 
occurs prior to detonation [17[. Since the great majority of the GW energy is radiated 
during the detonation itself, we expect both the DDT and DD models are likely to yield 
off-centered detonations whose gravitational wave signature is qualitatively similar to the 
GCD model prediction. 

Previous work explored the stochastic gravitational wave background due to the distribu- 
tion of unresolved "chirping" binary white dwarfs, and demonstrated that this background 
obscures the inflationary signal between frequencies of 10~^ Hz and 0.1 Hz 18|. The fre- 
quency range between 0.1 - 10 Hz was previously thought to be relatively free from grav- 
itational wave emission from astrophysical sources, with the primary contribution believed 



to originate from neutron star binaries 



19|, and possibly also core-collapse supernovae at 



cosmological distances, and Pop III (e.g., first-generation) stars [20|. In this paper, we con- 
sider an additional stochastic GW background source : SNe la at cosmological distances. 
If sufficiently strong, the SNe la stochastic background may be a source of noise to future 
detectors in the 0.1 - 10 Hz frequency range. Thus, the presence of a stochastic gravitational 
wave background due to SNe la may further impose additional challenges in extracting a 
possible inflationary GW background in this frequency range. 

The paper is organized as follows. In section II, we present the formalism underlying the 
computation of the cosmological gravitational wave background of stochastic sources. In 
section III, we discuss the cosmic rate of SNe la. In section IV, we review the gravitational 
wave signature of a single SNe la event predicted by the GCD mechanism. Finally, in 
sections V and VI, we present our results and conclude. 

STOCHASTIC BACKGROUND FORMALISM. 

The GW stochastic background is characterized by the dimensionless ratio of gravitational 
wave energy density to the critical energy density required to close the universe: 

(IcwU) = -P^,, (1) 

where pcw is the energy density of gravitational waves, / is frequency, and pc is the critical 
density, given as: 

Here G is the universal gravitational constant, and Hq is the Hubble parameter, which we 



take to be 70.4 km s ^ Mpc ^ 2l|. Using the expression for the energy density of supernovae 



GW sources, the gravitational wave energy density QcWif) can be expressed as an integral 



over redshift 



22|: 



fiaH.(/) = ^ r^-'%^% ^ (3) 



ScmiJo ~ E{z') df /.=/(!+.') 
where c is the speed of light, Rsn{z) is the SNe la event rate per comoving volume, dEg^/ df is 
the GW spectral energy density, and E{z) is a function of both redshift and the cosmological 
model : 

E{z) = ^nM{i + zf + nk{i + zf + nj,, (4) 



The spectral energy density is evaluated at the redshifted frequency fr = /(I + z). We take 
;he currently accepted cosmological parameters consistent with a A-cold-dark- matter model 



2l|: Qm = 0.27, Qa = 0.7, Qk = 0. 

The duty cycle is an important rate-dependent number that characterizes whether the 
stochastic background is discrete or continuous. It is determined by 

D= / dz'^V,{z)fGw{l + z'). (5) 



„ 

Here few is the typical rest-frame duration of the SNe la GW signal, which is set by the 
characteristic detonation time, and Vc{z) is the comoving volume as a function of redshift. 
If the duty cycle is less than unity, then the signal will consist of a series of rapid bursts 
separated in time, and relatively easily filtered from other signals, such as the primordial 
background of gravitational waves generated by inflation. In contrast, if the duty cycle is 
unity or larger, the background signal is continuous in time. Such a continuous signal acts 
as a foreground noise signal to a possible primordial background, and is more challenging 
to extract from the primordial background. 

THE COSMIC SNE lA RATE. 

The cosmic SNe la rate is observationally well-constrained only for low redshifts, z ^ 1. 
High-redshift rates are poorly-constrained due to their few-number statistics, but can be 



estimated using extrapolation and modeling. The SNe la rate is modeled as 23| : 

Rla[z{tc)] =11 r $ (tc - t;)p.[2;(t^)]< (6) 
Jtio 

here the star formation rate (SFR) per unit volume, is the progenitor-dependent delay 
time distribution (DTD), tz is the time corresponding to redshift z, is the current time, 
and 7] is the fraction of all stars which ultimately give rise to SNe la events. Mathematically, 
7] can be written as: 

7]= / Aia{M)aM)dM, (7) 

J Ml 

where ^(M) is the initial mass function of newly formed stars, Aia{M) is the SN la explosion 
efficiency as a function of initial stellar mass, and the integral is performed over the full range 
of initial stellar masses Mi — M2 of possible progenitors. If all SNe la progenitors exploded 
promptly at formation, the SNe la rate would simply follow the SFR. However because it 



takes time for a newly formed star to evolve to criticality, the rate is shifted from the SFR 
by the DTD. Since the SD and DD channels have different ranges of evolution times and 
relative frequencies that are uncertain, the modeled SNe la rate can diverge significantly 
due to the unknown DTD. 

As more high redshift SNe la data are collected, the DTD will be further contained 
and the progenitor model eventually determined. There are currently multiple methods for 
determining the DTD that use a range of techniques, environments, and redshifts (see 24 1 
for a review and comparison). Recent studies however generally point to a power-law DTD: 
$ oc with exponent /3 ^ — 1, suggesting a large component of DD events. Using new data 
collected from the Subaru Deep Field and previously published rates, j25| recently found 
(5 = —1.1 ± 0.1 (statistical) ±0.17 (systematic). 



THE TYPE lA GRAVITATIONAL WAVE SIGNAL. 

In our previous paper [l^ , we developed an estimate of the GW signal by post-processing 
3D simulations of the GCD mechanism within the SD channel of SNe la 26|. We found GW 



signal with characteristic frequency near 2 Hz, favoring the /i+ polarization. The strength 
of the GCD GW was found to be weak, but detectable as far as 1 Mpc by the Big Bang 
Observer (BBO), if the event happens to be oriented such that line of sight is close to the 
equator. In the current paper, we assume that all SNe la are generated through edge-lit 
detonations - including the GCD, DDT, and DD models - and that the GW spectral energy 
density (figure 4 of [14]) of the GCD model is typical of edge-detonated white dwarf (WDs) 
figure 4 [1^. For the current study, in order to accommodate high-redshift supernovae, we 
have extended the spectral density to higher frequencies. 

Summarizing our previous SNe la GW signal calculation, we simulated the hydrodynam- 
ical evolution of C/0 WD progenitors from the onset of defiagration through detonation 
in fully 3D numerical simulations [26]. The simulations were advanced within the FLASH 
framework with the Euler equations of inviscid, non-relativistic hydrodynamics, coupled to 
Poisson's equation for self-gravity, and an advection-diffusion reaction model of the thick- 
ened combustion front. FLASH is a modular, component-based application code framework 
created to simulate compressible, reactive astrophysical fiows and subsequently applied to 

These single-bubble GCD models produce intermediate 



a wide range of problems 



271, 



mass elements at a velocity coordinate ~ 11, 000 km /s, creating a layered structure of IME 
and Fe peak (NSE) products similar to observation 



Mm- However, single-bubble GCD 
models generally underproduce intermediate mass elements and overproduce ^^Ni. Conse- 
quently, the models are generally too luminous in comparison to Branch normal la events. 

Our simulation was initialized with a 1.37 Mq C/0 WD, with a single pre-ignited flame 
bubble of radius 16 km, and 40 km offset. As the WD evolved, we calculated, for each time 
step, the second-time derivative of the reduced quadrupole mass moments associated with 
the standard post-Newtonian expressions for the gravitational quadrupole radiation fleld 
from a 3D source 31|. 



RESULTS. 



To constrain the SNe la GW stochastic background, we consider an upper-bound SNe la 



rate associated with a prompt DTD model 23|, which could be characteristic of the shorter 



evolution of the SD channel, and the more probable /3 = —1.1 power-law DTD fit of [25 1. 
Here, we extend our integration to z = 10, which corresponds to the end of reionization, 
though our results do not depend sensitively upon this cutoff, as events at much higher 
redshift contribute negligible power to the overall background. 

We estimate the duty cycles associated with our chosen SNe la rates, assuming a charac- 
teristic detonation time of few ~ 0.5 s, obtained from the inverse of the characteristic 2 Hz 
detonation frequency [l^. Integrating to z = 10, we find duty cycles that are several orders 
of magnitude greater than unity for both assumed rates, so that the stochastic background 
is a continuous signal. 

To determine if the SNe la GW stochastic background will mask the inflationary sig- 
nal, we consider a scalar fleld (p inflationary model with slow-roll approximation. That is, 
the potential V{(j)) is assumed to be flat enough that the scalar fleld rolls slowly and the 
equations of motion governing inflation can be simplifled. We adopt the upper-bound esti- 
mate presented by [32], constrained by the energy scale of inflation V^^'^ set by the cosmic 
background radiation (CMB) and large scale structure (LSS) to be V'^^^ < 3.36 x 10^^ GeV. 

In flgure [T], we compare the sensitivities of third-generation detectors and the slow-roll 
inflation to our predicted GW stochastic background from cosmological GCD SNe la out 
to z = 10, assuming all are in the SD channel. Here we also present an upper-bound SNe 



la background estimate, derived by a simple estimate assuming that the SNe la event is 

n 

maximally asymmetric [14]. This upper-bound signal is obtained by shifting the prompt 
curve up by a factor of 100^, noting that flew cx; h'^ and the upper-bound amplitude is two 
orders of magnitude larger than our baseline calculated characteristic value. The baseline 
calculated value predicts that the stochastic background could be a source of noise for 
the inflationary signal for the ultimate form of the correlated Deci-Hertz Interferometer 



Gravitational wave Observatory (DECIGO) detector 
could be a source of noise for even the correlated BBO 
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3. 



The upper-bound estimate 



20|. 



CONCLUSIONS. 

A population of purely SD CCD SNe la at cosmological distances will generate a contin- 
uous GW stochastic background within the 1 - 2 Hz frequency range that could be a source 
of noise to third-generation detectors seeking the inflationary GW signal. The SNe la GW 
stochastic signal strength is relatively insensitive to the DTD, with the prompt DTD pro- 
ducing a slightly larger signal than the (3 = —1.1 power law. The SNe la background could 
be a significant source of noise for the ultimate DECIGO detector if the inflationary signal is 
at a level less than Qgw^o — 10^^^. However assuming a slow roll inflation that can be con- 
strained by the CMB/LSS, the inflationary signal could be as strong as Qcwho ^ 6 x 10^^^, 
in which case the SNe la background would be too weak to be a source of noise. If the SNe 
la GW signal strength is significantly stronger, such that the background is closer to the 
upper bound estimate, it could also be a source of noise to the correlated BBO. 

Previous work on the primordial inflationary gravitational wave spectrum has primarily 
bcused upon neutron-star binaries as a source of noise in the frequency range of 0.1 - 10 Hz 
lol . The gravitational wave signal generated by these binaries is relatively clean. Therefore 
it seems that it may be possible to subtract the neutron-star binary spectrum to reveal the 
underlying primordial gravitational wave spectrum. In contrast, the physics underlying the 
SNe la is far more complex, involving both hydrodynamics and nuclear combustion. Fur- 
thermore, though supernova rates are relatively well-constrained from optical observations 
at low redshift, they are poorly constrained at the high-redshifts responsible for the back- 
ground. In this respect, the numerous uncertainties surrounding the stochastic background 
of SNe la are similar to those of core-collapse SNe {2^. If the inflationary energy scale does 
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lie close to 10^^ GeV, the GW emission of SNe la will likely prove to be an additional source 
of noise to the primordial gravitational wave signal - one which will be very challenging to 
extract. 
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FIG. 1. Predicted and upper-bound SNe la GW stochastic background curves and the proposed 
correlated BBO and ultimate DECIGO detector sensitivities scaled by the dimensionless Hubble 
factor /iQ = 0.72. Two possible slow-roll inflationary GW stochastic background curves are pre- 
sented, an upper-bound estimate derived from the CMB and a second estimate close to the level 
of the prompt SNe la background. 
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